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The renormalized trajectory (RT) is determined from two different Monte Carlo renormalization group tech- 
niques with (5-function block spin transformation in the multi-dimensional coupling parameter space of the two- 
dimensional non-linear sigma model with 0(3) symmetry. At a correlation length ^ 3-5, the RT is shown to 
break away from the straight line of the fixed point trajectory (FPT) which is orthogonal to the critical surface 
and originates from the ultraviolet fixed point (UVFP). The large N calculation of the RT is also presented in 
the coupling parameter space of the most general bilinear Hamiltonian. The RT in the large N approximation 
exhibits a similar shape with the sharp break occurring at a somewhat smaller correlation length. 



1. INTRODUCTION 

A sketch of the RT in the phase diagram of 
asymptotically free (AF) field theories is shown in 
Fig. 1 as plotted in the infinite dimensional space 
of coupHngs Ki. The inverse coupHng 1/Ki of 
the standard lattice action is singled out to label 
the horizontal axis. The critical manifold is given 
in Fig. 1 by the Ki = oo plane and UVFP desig- 
nates the AF fixed point on the critical surface. 
The continuum theory at finite lattice correlation 
length is defined by the RT which fiows along the 
unstable direction from the UVFP to the high 
temperature fixed point (HTFP) with vanishing 
correlation length at Ki = [1]. 

It has been known for a long time that, in con- 
trast with the standard lattice action, each point 
on the RT defines a perfect lattice action [2] which 
is free of cutoff effects at any finite correlation 
length [1] . The determination of this perfect lat- 
tice action, however, has not been considered fea- 
sible in the past, mainly because a reahstic ap- 
proximation to the RT in the infinite dimensional 
coupling parameter space was thought to require 
the uncontrolled truncation of a large number of 
terms in the blocked Hamiltonian. 

Considerable progress has been made recently 
by Hasenfratz and Niedermayer [2] who realized 
that the fixed point lattice action of the UVFP 
can be determined from a classical saddle point 
problem in AF field theories. The FPT which 
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originates from the UVFP is suggested to be a 
good approximation to the RT at sufficiently large 
correlation lengths [2] (dashed line at if i < oo in 
Fig. 1.) It was also demonstrated that the fixed 
point action (FPA) can be rendered very short- 
ranged by optimizing the block-spin renormaliza- 
tion group transformation with the expectation 
that an approximate FPT which is Vjased on a 
truncated FPA, with the very small long range 
couplings neglected, will exhibit almost no cut- 
off dependence at large correlation lengths. 

It remains an important question to deter- 
mine tlie crossover region in the correlation length 
where the RT departs from the FPT. The per- 
formance of the approximate FPT has first been 
tested by Hasenfratz and Niedermayer in the 2d 
non- linear 0(3) sigma model which is known to 
be AF [2]. They carried out a pilot study using a 
truncated FPT with 24 difl^erent couplings. Sev- 
eral tests showed that the residual cutoff effects 
were not visible even down to a correlation length 
of five suggesting that truncation effects are neg- 
ligible. This indirectly implies that the FPT runs 
close to the RT in an extended range of the lattice 
correlation length. 

The RT eventually has to break away from the 
FPT since it flows into the HTFP as depicted in 
Fig. 1. For practical applications, it is important 
to determine the crossover region in the correla- 
tion length where the RT breaks away from the 
FPT. In this contribution we determine numeri- 
cally the position of the RT in a flnite dimensional 
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Figure 1. Phase digram of A F field theories. 

subspace of the infinite dimensional coupling pa- 
rameter space and compare it with the FPT. 

2. THE FIXED POINT ACTION 

The path integral of the 0(N) non-linear 
sigma model on the lattice is given by 
Z = JY\^'D(f)xexp{—/3'H), where x is a two- 
dimensional lattice vector to label the sites, (j)x 
is an A'' component scalar field, and J'D(j)x = 

f-^ Tlx d(j>xS{l — (j)x) is the 0(N) invariant mea- 
sure. The lattice hamiltonian Ti is not unique. 
The only constraints come from the symmetries 
of the system and the fact that it ought to reduce 
in the classical continuum limit to the continuum 
hamiltonian. The lattice hamiltonian with these 
constraints can be parametrized as 

n=-H2+n4 + -H6 + ... , (1) 

r X 

"^4 =^Cx-,,X2,X3,X4 (1 - (I^Xl<i>X2)i'^ " '^^Xs^'xi), (3) 

where 7^2 includes all terms that are bilinear and 
'H4 , . . . designate all the other terms which arc 
quartic or of higher order in the field variables. 
The summations in Eqs. (2) and (3) are over lat- 
tice vectors r, x, xi, x-2. x?,, and x^. The standard 
lattice action is given by the first two terms in 
Eq. (2) with r = (1, 0), (0, 1), and = /3p(ri). 
Let us now consider the following J-function 



Table 1 

The six largest couplings p*{r). 



r 


p*{r) 


ri 


= (1,0) 


-3.42839 


r2 


= (2,0) 


-fO. 74981 


rs 


= (1,1) 


-1-0.32486 




= (3,0) 


-0.16871 




= (4,0) 


+0.03877 


re 


= (3,1) 


-0.01976 



block spin transformation, 

where a blocked lattice site x' is assigned to a 
2x2 cell of sites x on the unblocked lattice. For 
details on how the FPA can be computed in AF 
field theories we refer the reader to Ref. [2]. The 
authors of Ref. [2] showed that the fixed point 
couplings p* (r) agree with the ones obtained ear- 
lier for the non-interacting model with Gaussian 
block spin transformation [3], 



(6) 



P*(p) 



+ 00 



2 



{pi + 27rni)2 -I- (p2 + 27rn2)2 



n 



sin"(p»/2) 

\ (pi/2 + miT) 



(7) 



which holds independently of N . We have deter- 
mined the coupHngs for several lattice vectors r 
by evaluating Eqs. (6) and (7) numerically. The 
values of the six largest couplings are given in 
Table 1. The interactions are indeed very short- 
ranged since the couplings decrease rapidly when 
the distance \r\ between two spins grows. 

3. LARGE N CALCULATION 

In the large N limit, Hirsch and Shenker de- 
rived a recursion relation for the blocked two- 



3 



-0.14 




-0.05 - 



ln(l/f) 

Figure 2. Flow lines of p{r2) I p{r\) in the large 
N limit. 

point function in momentum space [4], 



G'(P') 



^Ep',G{f +7r0 +7r0 ' 



^(^^') = n 



sin^ p' 



(8) 
(9) 



where I is a vector with components equal to 
or 1. designates the number of points on the 
lattice before blocking, and G{p) = [I3p{p) + X]~^ 
is the saddle point solution with the mass gap \/X 
determined from the gap equation, -p- X]p ^(p) = 
1. Eq. (8) can now be iterated to determine 
the blocked propagators after repeated block- 
spin transformations. In the subspace of bi- 
linear Hamiltonians we can determine the cou- 
plings (3p{r), (3'p'{r), etc. from the inverse 
propagators by Fourier transformation (3p{r) = 



etc. 



We have iterated Eq. (8) numerically, using two 
different expressions for the unblocked propaga- 
tor. The first choice originates from the bilinear 
Hamiltonian that includes the four largest cou- 
plings in Table 1 and the second choice is the 
fixed point propagator resulting from Eq. (7). In 
the first case we start the fiow lines, which after 
several blocking steps approach the RT, at dif- 
ferent positions in the four-dimensional coupling 
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Figure 3. Flow lines of p{r2) I p{r\) for N 



parameter space. In the second case, the flow 
lines start on the FPT and it is interesting to see 
how the blocked couplings break away from the 
FPT after a few block spin transformations. 

As an example, we have displayed in Fig. 2 the 
flow lines of the ratio p{r2) / p{ri) as a function of 
ln(l/^), with ^ designating the correlation length. 
The plots for the other couplings look very sim- 
ilar and we refer the reader for more details to 
Ref. [5]. The UVFP is located in this plot at 
ln(l/0 = -00 and the HTFP at ln(l/0 = +oo. 
The ratio p* {r2) / p* (ri) has been represented in 
this graph by the horizontal dashed line. The flow 
fines which start from the four-coupfing Hamil- 
tonian are represented by crosses. Fig. 2 shows 
that those flow lines that start at large correlation 
lengths approach the FPT within a few block spin 
steps indicating that the FPT is indeed a very 
good approximation to the RT at large correla- 
tion lengths. However, a different behavior is seen 
at small correlation lengths where the flow fines 
approach the RT which runs into the HTFP with 
the ratio p{r2) / p{ri) vanishing, as can be shown 
by the high temperature expansion [5]. Fig. 2 
shows that the sharp break from the FPT occurs 
at a correlation length of 2 — 3. The three fiow 
lines that start on the FPT are represented in the 
graph by squares, full, and open circles. 
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4. NUMERICAL SIMULATION OF THE 
0(3) MODEL 

The renormalization group flows in the pre- 
vious section were restricted to the subspace of 
bilinear Hamiltonians in the large N limit. In 
contrast, the Monte Carlo renormalization group 
method will allow us to study the exact flows of 
selected couplings in the 0(3) model. Although 
interaction terms which do not fit on the finite 
lattice arc truncated, their effects are expected to 
be negligible. The only practical limitation on the 
method is the signal to noise ratio in measuring 
very smafl blocked coupHngs. 

The 1-clustcr algorithm has been used to gen- 
erate unblocked spin configurations on a 256 x 256 
lattice using a Hamiltonian with the first four 
coupHngs of Table 1. The spin configurations 
have been blocked four times down to a lattice 
of size 16 X 16 using the block spin transforma- 
tion of Eq.(4). Microcanonical [6] and canonical 
[7] demon methods have been implemented to in- 
fer the blocked couplings from the blocked spin 
configurations. For more details and a compar- 
ison of the two methods, we refer the reader to 
ref. [5]. 

To illustrate our results, we have plotted again 
in Fig. 3 the renormalization group fiow of the 
ratio p{r2)/p{ri) as a function of ln(l/^). Most 
of the data have been obtained with the micro- 
canonical method (crosses). We have also in- 
cluded two renormalization group fiow lines which 
have been generated by the canonical demon 
method (circles), one of which with the same ini- 
tial couplings as a selected microcanonical fiow 
line. Fig. 3 shows that the two fiow fines coin- 
cide and the fiow pattern is very similar to the 
one seen in Fig. 2. The only difference is that the 
sharp break from the FPT occurs now at a some- 
what larger correlation length. The RT of the 
0(3) model is traced by the grey line. For com- 
parison, we have also included in the graph the 
RT in the large A'' fimit (heavy fine) (cf. Fig. 2). 

5. CONCLUSIONS 

We have demonstrated in this paper by explic- 
itly computing the RT in the 0(3) model that 



the FPT provides a good approximation to the 
RT in the region of the coupfing parameter space 
where the correlation length is large. A significant 
break occurs, however, at a correlation length of 
about 3-5 where the RT sharply departs from the 
FPT and fiows into the HTFP. The RT in the 
large N fimit exhibits a similar behavior with the 
sharp break shifted to somewhat smafier correla- 
tion length. Although wc have only shown the 
projection of the exact RT in the K1-K2 plane, 
several other projections of the fiows were deter- 
mined as well [5]. 

Both ends of the RT can be determined with- 
out the utilization of the Monte Carlo renormal- 
ization group. The RT in the large correlation 
length regime, near the UVFP, is well approxi- 
mated by the FPT, whfie the RT in the regime 
of small correlation lengths, near the HTFP, can 
be obtained from a high temperature expansion. 
It is likely that both approximations will break 
down in the crossover region where the Monte 
Carlo renormalization group technique may re- 
main the only useful tool. It is not clear whether 
this regime wifi be important in practical appli- 
cations. 

The studies we presented will be useful to ex- 
tend to non-Abelian lattice gauge theories in four 
dimensions [8]. 
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